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Abstract 

Robust heteroclinic cycles in equivariant dynamical systems in have been a sub¬ 
ject of intense scientific investigation because, unlike heteroclinic cycles in they can 
have an intricate geometric structure and complex asymptotic stability properties that 
are not yet completely understood. In a recent work, we have compiled an exhaustive 
list of finite subgroups of 0(4) admitting the so-called simple heteroclinic cycles, and 
have identified a new class which we have called pseudo-simple heteroclinic cycles. By 
contrast with simple heteroclinic cycles, a pseudo-simple one has at least one equilib¬ 
rium with an unstable manifold which has dimension 2 due to a symmetry. Here, we 
analyse the dynamics of nearby trajectories and asymptotic stability of pseudo-simple 
heteroclinic cycles in 


1 Introduction 

It is known since the 80’s that vector fields defined on a vector space V (or more generally a 
Riemannian manifold), which commute with the action of a group F of isometries of V, can 
possess invariant sets which are structurally stable (within their F symmetry class) and which 
are composed of a sequence of saddle equilibria ... ,^m and a sequence of trajectories Ki, 
such that Kj belong to the unstable manifold of as well as to the stable manifold of .^,+1 
for all z = where ^m+i = G- These objects are called robust heteroclinic cycles 


1 


and under certain conditions they can be asymptotically attracting. Robust heteroclinic 
cycles have been considerably studied in low-dimensional vector spaces ms]. Although 
their properties and classihcation are well established in R = the case of appears 
to be much richer and yet not completely investigated. In | 8 ] the so-called simple robust 
heteroclinic cycles were introduced, which can be dehned as follows. Let Aj denote the 
isotropy subgroup of an equilibrium ^j. The heteroclinic cycle is simple if (i) the hxed point 
subspace of Aj, denoted Fix (Aj), is an axis; (ii) for each j, is a saddle and is a sink in 
an invariant plane Fix(Sj) (hence Sj C Aj flAj+i); (hi) the isotypic decompositioi]0 for the 
action of Aj in the tangent space at only contains one-dimensional components. In [13] 
we have found all hnite groups F cO(4) admitting simple heteroclinic cycles (i.e. for which 
there exists an open set of F-equivariant vector fields possessing such an invariant set). We 
have also pointed out that the dehnition of a simple heteroclinic cycle in former works m 
had omitted condition (iii), which was implicitly assumed. If this condition is not satisfied, 
then the behaviour of the heteroclinic cycle turns out to be more complex, because there 
exists at least one equilibrium at which the unstable manifold is two dimensional and, 
moreover, is invariant under a faithful action of the dihedral group for some k > 3. We 
called heteroclinic cycles satisfying (i) and (ii) but not (iii) pseudo-simple. 

The aim of the present work is to investigate the dynamical properties of pseudo-simple 
heteroclinic cycles in We shall therefore consider equations of the form 

X = /(x), where /(yx) = 7/(x) for all 7 G F, F C 0(4) is hnite (1) 

and / is a smooth map in and which possess a pseudo-simple cycle. 

Our main result stated in Theorem 1, Section |31 is that if F C S0(4), then the pseudo¬ 
simple cycle is completely unstable. Namely, there exists a neighborhood of the heteroclinic 
cycle such that any solution with initial condition in this neighborhood, except in a subset 
of zero measure, leaves it in a hnite time m- Then we aim at analysing more deeply what 
the asymptotic dynamics can be in a neighbourhood of the heteroclinic cycle. We do this 
by focusing on specihc examples. In Section 0] we introduce a subgroup of S0(4) which is 
algebraically elementary, however possessing enough structure to allow for the existence of 
pseudo-simple heteroclinic cycles. This group is a (reducible) representation of the dihedral 
group D3. Hence, F is algebraically isomorphic to D3. In this example, the cycle involves 
two D 3 invariant equilibria, .^1 and ^ 2 , such that for both the linearization df{^i) possess 
a double eigenvalue with a D 3 invariant associated eigenspace, one being negative and the 
other being positive (unstable). We show that when the unstable double eigenvalue is small, 
an attracting periodic orbit can exist in the vicinity of the F-orbit of the heteroclinic cycle, 
the distance between the two invariant sets vanishes as the unstable eigenvalue tends to 
zero. This result is numerically illustrated by building an explicit third order system for 
which a pseudo-simple cycle exists. In the next Section 0] we extend the group F by adding 
a generator which is a reflection in The resulting group is isomorphic to D 3 x Z 2 and 
possesses the same properties for the existence of a pseudo-simple cycle. We show that in 

^The isotypic decomposition of the representation of a group is the (unique) decomposition in a sum of 
equivalence classes of irreducible representations. 
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this case the cycle is not completely unstable anymore, but asymptotically fragmentarily 
stable im instead. This means that the cycle is not asymptotically stable, however in its 
any small neighbourhood there exists a set of positive measure, such that any solution with 
initial condition in that set remains in this neighbourhood for alH > 0 and converges to the 
cycle as f —)■ oo. This result is then illustrated numerically. 

In Section [ 6 ] we present another example of a group in SO(4) admitting pseudo-simple cycles 
with a more complex group structure and acting irreducibly on The Theorem of existence 
of nearby stable periodic orbits proven in Section H] holds true for the same reasons as in the 
D 3 case. We then investigate numerically the asymptotic behavior of trajectories near the 
cycle and show that periodic orbits of various kinds can be observed as well as more complex 
dynamics when the unstable double eigenvalue becomes larger. 

In the last Section of this paper we discuss our hndings and possible continuation of the 
study. 

2 Pseudo-simple heteroclinic cycles 

In this Section we introduce the notion of a pseudo-simple heteroclinic cycle and notations 
which will be used throughout the paper. 

2.1 Basic definitions 

Given a hnite group T of isometries of let S be an isotropy subgroup of T (the subgroup 
of elements in T which fix a point in M'^). We denote by Fix (E) the subspace comprised of 
all points in which are fixed by S. 

Note that (i) if S C A C T, then Fix (A) C Fix(E); (ii) 7 ■ Fix (E) = Fix( 7 E 7 "^) for any 
7 G F. 

In the following we shall assume the F-equivariant system ([T]) admits a sequence of isotropy 
subgroups Aj C F, j = 1,..., M, such that the following holds: 

(i) dimFix(Aj) = 1 (axis of symmetry). We denote Lj = Fix(Aj). 

(ii) On each Lj, there exists an equilibrium of ([T]). 

(ii) For each j, Lj = Pj-i O Pj, where Pj is a plane of symmetry: Pj = Fix (Ej) for some 
Ej C Aj. We set Pm = Pq- 

(iv) is a sink in Pj-i and a saddle in Pj. Moreover, a saddle-sink heteroclinic trajectory 
Kj of ([T]) connects to in Pj (j = 1 ,..., M). Note that km connects to in 
Pm- 

Conditions (i) - (iv) insure the existence of a robust heteroclinic cycle for ([T]). We denote 
by Jj the Jacobian matrix df{^j). Since Lj is flow-invariant, Jj has a radial eigenvalue rj 
with eigenvector along Lj and by (iv) we have rj < 0. Moreover (iv) implies that in Pj-i, Jj 
has a contracting eigenvalue —Cj {cj > 0 ), while in Pj it has an expanding eigenvalue Cj > 0 
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corresponding to the unstable eigendirection in that plane. The remaining eigenvalue tj is 
called transverse, its associated eigenspace is the complement to © Pj in 

Definition 1 We say that the group T admits robust heteroclinic cycles if there exists an 
open subset of the set of smooth T-equivariant vector fields in R", such that vector fields in 
this subset possess a (robust) heteroclinic cycle. 

In order to insure the existence of a robust heteroclinic cycle it is enough to hud m < M 
and 7 G r such that a minimal sequence of robust heteroclinic connections ^ ^ ^m+i 

exists with ^m+i = 7 ^i (minimal in the sense that no equilibria in this sequence satisfy 
for some 7 ' G T and 1 < i,j < m). It follows that 7 ^ = 1 where /c is a divisor of 

M. 

Definition 2 The sequence ^ • • • —t and the element 7 define a building block of the 

heteroclinic cycle. 

The asymptotic dynamics in a neighborhood of a robust heteroclinic cycle is what makes 
these objects special. A heteroclinic cycle is called asymptotically stable if it attracts all 
trajectories in its small neighborhood. When a heteroclinic cycle is not asymptotically 
stable, it can however possess residual types of stability. Below we dehne the notions, which 
are relevant in the context of this paper, see HU. Given a heteroclinic cycle X and writing 
the flow of ([T]) as the 5-basin of attraction of X is the set 

Bs{X) = {x G R^ ; d{^t{x), X) < 6 for alH > 0 and lim (i(<hi(a:), X) = 0}. 

t—^+00 


Definition 3 A heteroclinic cycle X is completely unstable if there exists 5 > 0 such that 
Bs{X) has Lebesque measure 0. 

Definition 4 A heteroclinic cycle X is fragmentarily asymptotically stable if for any 5 > 0, 
Bs{X) has positive Lebesque measure. 

2.2 Simple and pseudo-simple heteroclinic cycles 

Recall that the isotypic decomposition of a representation T of a (hnite) group G in a vector 
space V is the decomposition V = ® • • • ® where r is the number of equivalence 

classes of irreducible representations of G in R and each = Tjy, is the sum of the 
equivalent irreducible representations in the j-th class. This decomposition is unique. The 
subspaces are mutually orthogonal (if G acts orthogonally). 

Then the following holds [T3] . 

Lemma 1 Let a robust heteroclinic cycle in R^ be such that for all j: (i) dimP^ = 2, (ii) 
each connected component of Lj \ {0} is intersected at most at one point by the heteroclinic 
cycle. Then the isotypic decomposition of the representation of Aj in R^ is of one of the 
following types: 
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1. Lj ©-‘- Vj ©■‘“ Wj ©■‘“ Tj (the symbol ©■‘“ indicates the orthogonal direct sum). 

2. Lj ©-‘- Vj ©-‘- Wj where Wj = Wj © Tj has dimension 2. 

3. Lj ©■‘“ Vj ©■*■ Wj where Vj = Vj ® Tj has dimension 2. 

In cases 2 and 3, acts in Wj (respectively, Vj) as a dihedral group in for some 
m > 3. It follows that in case 2, Cj is double (and Cj = tj) while in case 3, —Cj is double 
(and —Cj = tj). 

Definition 5 A robust heteroclinic cycle in satisfying the conditions 1 and 2 of Lemma 
[H is called simple if case 1 holds true for all j, and pseudo-simple otherwise. 

Remark 1 DefinitionlE implies that any pseudo-simple heteroclinic cycle in a T-eguivariant 
system m, where V C SO(4), has at least one I < j < M such that Ej = Zk with k > 3. 

Remark 2 An order two element a in SO(4) whose fixed point subspace is a plane P must 
act as —Id in the plane fully perpendicular to P. Nevertheless, to distinguish it from 
other rotations fixing the points on P, we call a a plane reflection. 

3 Instability of pseudo-simple cycles with P C 50(4) 

In this Section we prove the following Theorem: 

Theorem 1 Let be a pseudo-simple heteroclinic cycle in a T-eguivariant 

system, where T C S0(4)- Then generically the cycle is completely unstable. 

Proof: Following |H1 [9l EH |I2], to study asymptotic stability of a heteroclinic cycle, we 
approximate a ’’first return map” on a transverse (Poincare) section of the cycle and consider 
its iterates for trajectories that stay in a small neighbourhood of the cycle. 

In Section 12.11 we have defined radial, contracting, expanding and transverse eigenvalues of 
the linearisation Jj = df{f,j). Let {u,v,w,q) be coordinates in the coordinate system with 
the origin at fj and the basis comprised of the associated eigenvectors in the following order: 
radial, contracting, expanding and transverse. If S is small, in a ^-neighbourhood of fj the 
system (E]) can be approximated by the linear system 

ii = —rjU 
V = —CjV 
W = CjW 

q = tjq. 

Here, {u,v,w,q) denote the scaled coordinates {u,v,w,q) = {u,v,w,q)/S. In fact a version 
of the Hartman-Grobman theorem exists |2], which allows to linearize the system with a 
{k > 1) equivariant change of variables in a neighbourhod of fj if conditions of nonresonance 
are satisfied between the eigenvalues (which is a generic condition here). 


( 2 ) 
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Let (uo,Vo) be the point in -Pj-i where the trajectory kj^i intersects with the circle 
= 1 , and h be local coordinate in the line tangent to the circle at the point 
We consider two crossections: 

= {(h,w,q) : \h\,\w\,\q\<l} 

and ^ 

= {(n,n,'«;,g) : |m|, |n|, |g| < 1, w = 1}. 

Near trajectories of (II]) with initial condition in hit the outgoing section 
after a time r, which tends to inhnity as the initial condition comes closer to Pj-i- The 
global map 'ijjj : —>■ associates a point where a trajectory crosses with the 

point where it previously crossed It is dehned in a neighborhood of 0 (in the local 

coordinates in and is a diffeomorphism. The hrst return map is now dehned as the 

composition 'g = qm ° - ^ gi, where Qj = xfjj o 0^-. 

As it is shown in mm, in the study of stability only the coordinates that are in 
are of importance. By (pj and tpj we denote the maps pj and tpj restricted to P^ (to be more 
precise = Pj-_i n and = P/ n 

Now, by dehnition a pseudo-simple cycle has at least one connection Kj : such 

that Sj = Zfc with k > 3 (see remark [1]). We can assume that j = 1. The equilibria and 
^2 belong to the plane Pi = Fix (Ei). We now prove existence of £ > 0 such that 

02 o o 0i(p|*''^(e:)) n = 0 , (3) 

where 

= {{w,q) E : \{w, q)\ < e} and = {{v, q) E : \{v,q)\<e}. 

Evidently, because of (I3|) the hrst return map in cannot be completed in general, 

which implies that the cycle is completely unstable. 

To prove (I3|) we employ polar coordinates (pi, 6 *i) and (p 2 , 6 * 2 ), in and respec¬ 

tively, such that Vi = pi cos6*i, qi = pi sin6*i, W 2 = P 2 cos 02 and q 2 = P 2 sin 02 - Due to (I2|) in 
the leading order the maps 0i and 02 are 

(pi,0i) = 0i(wi,gi) = (no,iwJ'^''\arctan(gi/no,i)), (4) 

(^2,02) = 02(P2,02) = (no,2(p2COS02)^=/®^tan02). (5) 

In the leading order the global map 0i : -E- D linear. The map commutes 

with the group S 2 acting on P 2 as rotation by 27r/fc, where k > 3, therefore in the leading 
order 0i is a rotation by a hnite angle composed with a linear transformation of p, 

(p2, 02) = 'Piipi, 0i) = {Api, 01 + 0), (6) 

where generically 0 7 ^ Nn/k. 
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Denote 


a= min IQ — Nn/kl 

l>N>2k 


(7) 


and set 


0 < £ < min 



a 

no,i tan - 


( 8 ) 


Any {wi,qi) G satisfies qi < e, therefore dl]) and ([HD imply that 6i < a/2. Hence, 

due to dHD-dSD 16^2 — Nn/k\ > a/2 for any k. The steady state ^2 has k symmetric copies 
(under the action of symmetries cr G E 2 ) of the heteroclinic connection K 2 ■ ^2 ^ ^3 which 
belong to the hyperplanes 6*2 = Nn/k with some integer N^s. Due to (jSD and (jHD, the distance 
of {v 2 ,q 2 ) to any of these hyperplanes is larger than tan(a/2), which implies (jHD. QED 


4 Existence of nearby periodic orbits when P C /S'0(4): 
an example 

We have proved that any pseudo-simple heteroclinic cycle in a T-equivariant system, where 
r cS0(4), is completely unstable. In this Section we show that nevertheless trajectories 
staying in a small neighbourhood of a pseudo-simple cycle for all f > 0 can exist. Namely, 
we give an example of a subgroup of SO(4) admitting pseudo-simple heteroclinic cycles and 
prove generic existence of an asymptotically stable periodic orbit near a heteroclinic cycle if 
an unstable double eigenvalue is sufficiently small. 

4.1 Definition of P and existence of a pseudo-simple cycle 

We write (xi,|/i,X 2 , 2 / 2 ) G and Zj = Xj - 1 - iyj. Let T be the group generated by the 
transformations 

P ■ {zi,Z2) H- {zi,e^Z2), k{zi,Z2) H- (zi,Z2) (9) 

This group action obviously decomposes into the direct sum of three irreducible representa¬ 
tions of the dihedral group Ds: 

(i) the trivial representation acting on the component xi, 

(ii) the one-dimensional representation acting on yi by Kyi = —yi, 

(iii) the two-dimensional natural representation of D 3 acting on Z 2 = {x 2 ,y 2 )- 

There are two types of hxed-point subspaces and one type of invariant axis for this action: 

(i) Pi = {(xi, yi, 0, 0)} = Fix (p), 

(ii) P 2 = {(xi,0,X2,0)} = Fix(K), 

(iii) L = Pi n P 2 = Fix (D 3 ). 

Note that Pi is F-invariant, while P 2 has two symmetric copies, P 2 = pP 2 and P!/ = p^P 2 - 

Proposition 1 The group F, generated by p and k admits pseudo-simple cycles (see 
Definition\^. 
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Proof: The proof of existence of an open set of smooth G eqnivariant vector helds with 
saddle-sink orbits in Pi and in P2 connecting eqnilibria ^1 and ^2 lying on the axis L and snch 
that ^1 is a sink in P2 and ^2 is a sink in Pi, goes along the same argnments as in Lemma 
5 of [l 3 ]. We just need to check that the cycle is pseudo-simple. This comes from the fact 
that the Jacobian matrix of the vector held taken at .^1 or ,^2 has a double eigenvalue with 
eigenspace {(0, 0, X2,1/2)} due to the action of p, which generates the subgroup Z3 (negative 
eigenvalue at .^1 and positive eigenvalue at .^2)- QED 

By letting T act on the heteroclinic cycle, one obtains a 6-element orbit of cycles which 
all have one heteroclinic connection in Pi and the other one in P2, P2 or P^'. 

4.2 Existence and stability of periodic orbits 

We prove the following Theorem. 

Theorem 2 Consider the T -equivariant system 

X = /(x,/i), where /(7x,/i) = 7/(x,/i) for all 7 G T, (10) 

/ : X R —)■ is a smooth map and T is generated by p and k Let ^1 and ^2 be the 
equilibria introduced in proposition{l\ —cj and ej be the eigenvalues of df{fj), such that —ci 
and 62 are double eigenvalues with a natural action of D3 in their eigenspaces and ei and 
—C2 are the eigenvalues associated with the one-dimensional eigenspace, where the action of 
T in not trivial. Suppose that there exists po > 0 such that 

(i) 62 < 0 for —/io < p < 0 and 62 > 0 for 0 < p < po; 

(a) for any 0 < p < po there exist heteroclinic connections Ki = (hW(^i)nPi) 71114(^2) 7^ 0 
and K2 = (Wu{^2) H P2) fl ITs(^i) 7^ 0 , introduced in proposition{l\ 

Denote by X the group orbit of heteroclinic connections Ki and K2: 

X = (U.y6r7^i) U(U7er7^2)- 


Then 

(a) if 3ci < Cl then there exist p' > 0 and J > 0, such that for any Q < p < p' almost all 
trajectories escape from Bs{X) as t —)■ 00; 

(b) if 3 ci > Cl then generically there exists periodic orbit bifurcating from X at p = 0 . To be 
more precise, for any 6 > 0 we can find p{6) > 0 such that for all 0 < p < p{6) the system 
/ fldl) possesses an asymptotically stable periodic orbit that belongs to Bs{X). 

The rest of this Section is devoted to the proof of the Theorem. We start with two 
Lemmas which describe some properties of trajectories of a generic D3-equivariant systems. 


4 . 2.1 Two Lemmas 


A generic Da-equivariant second order dynamical system in C is 

z = az^ f3z'^. 

In polar coordinates, z = re'®, it takes the form 

r = ar + ( 3 r‘^ cos 39 , 

9 = —/Sr sin 36 *. 


( 11 ) 


( 12 ) 


We assnme that a > 0 and /S > 0 . The system has three invariant axes with 9 = Ktt/S, 
= 0 , 1 , 2 , and three eqnilibria off origin, with r = a/f 3 and 9 = { 2 k + l)(/r/ 3 , k = 0 , 1 , 2 . 
We consider the system in the sector 0 < 6* < 7r/3, the complement part of C is related to 
this sector by symmetries of the gronp ^3. 

Trajectories of the system satisfy 


dr 


a + /Sr cos 39 


/S sin 39 


(13) 


Re-writing this eqnation as 


dr cos 39 

d6' sin 36 * 


a 


/S sin 36 * ’ 

mnltiplying it by sin^'^^ 36 * and integrating, we obtain that 


r sin^/^ 39 


rO 

--S{ 9 ) + C, where ^( 0 ) = / sin’^/^ 3 M 0 , 

P Jo 


( 14 ) 


which implies that 

r sin^/^ 36 * -h ^S{ 9 ) = sin^/^ 36 *o -6 ^*S'(6*o) ( 15 ) 

/J fJ 

for the trajectory throngh the point (ro,6*o). Let 6*(r;ro,6*o) be a solution to flT^ with the 
initial condition r( 0 ) = tq and 9 ( 0 ) = 6*0. Then 


C = ^S{ 9 ), where 6* = lim 6*(r; rg, 6*0). 

P r^O 


We can re-write flTT)) as 



C\ c\ ~ 

r sin^/^ 36* = — — S{ 9 ) + 

/J jJ 

(16) 

Note, that 


S( 5 )^ 

rfi/6) 

-;—, where T is the Gamma function. 

3 r(2/3)’ 

(17) 


Lemma 2 Given j 3 > 0 and £ > 0 , there exists cto > 0 such that for any 0 < a < ag, 
0 < vq < e and 0 < 6*0 < 7r/3 


r sin 6*(r; ro, 6*0) < e for all t > 0. 


(18) 


9 







0 


Figure 1: Trajectories of the system ffT^ in the sector 0 < 6 * < vr/S. 

Proof: Set oq = £/^/4. From ffl^ 

—rsin^ = r sin 6 ' + Or cos 9 = rsin^fa — (3r cos9/2). (19) 

at 

From ffTTj) . for 0 < 6^0 < vr/3 we have lim^^oo 0{r; ro, Oq) = 0, which implies that d/dt(r sin 6 *) < 
0 for large r. Therefore, the maximum of rsin 6 ' is achieved either at t = 0, or at the point 
where 2a = (3rcos6. Since 0 < 0 < vr/S, at this point rsin 6 * = 2asm.6/{(3 cosO) < Aa/j3 < e. 
Since tq < e, we have r sin 6 < e att = 0 as well. (The behaviour of trajectories of the system 
dUis shown in Fig. [T]) QED 


Lemma 3 Let T{ro, Oq) denotes the time it takes the trajectory of the system starting 
at (ro, 6^0) to reach r = 1 and 'd(ro, 6*0) denotes the value of 9 at r = 1. Then 


(i) r(ro,0) satisfies 

gar(ro,0) _ ^0 + 

ro(l + a/(3)' 


(a) r(ro,6'o) satisfies 


r(ro, 9o) > r(ro, 0) for any 0 <9o<7i/3. 


( 20 ) 


(Hi) 'd(ro,6*o) satisfies 

, (y ! Ol 

sin^/^3'd(ro,6'o) +-5'(79(ro,6'o)) = rosin^^^ 39o + -S{9o). (21) 

(iv) Given C > 0, /3 > 0 and 0 < Oq < n/S, for sufficiently small a and uq 

<^9(ro,0o). 
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(v) For small a, tq a, r' tq and 6*' -C 1 


|g-Cr(ro+r',0o+0') _ ^-CTiroMj < g-Cr(ro.0o) | ^ + 30'^(0o)r (rQ, ^o) | • 

aro 


Proof: To obtain (i), we integrate the first eqnation in ffT^ . where we set 0 = 0. Since in 
m» < 6' implies r{r,6) > r{r,d'), (ii) holds trne. The eqnality (hi) follows from flT5|l . 

Since 9 < 0 (see flT^ ), (iii) implies that '9{ro,9o) -C 1 for small a and tq. Therefore 
^9(?^o,6'o) ~ ^(^’^osin^/^36'o + ^S'(6»o)^ > ^ ^max(a,/3ro) min(sin^/^ 300, S'(6»o))^ . (22) 
Due to (ii), to prove (iv) it is sufficient to show that 

e-^^('-O’O) <^(ro,0o). 

Because of (i), for small ro and a 


-Cr(ro,0) _ f ro(l + a/^) \ / ~ [(1 _ g)l/s]C/(a+/3ro) ^ _ 

\ ro +a/P J ^ ^ 


a 


ro + a/(3 

Since supo<s<i(l “ we have 


a + f3ro 


g-CT(ro,0) ^ ^-C/ max(a,/3ro) 


(23) 


Together, fl2^ and (123|1 imply (iv). 

For r ^ a the first equation in (IT^ can be approximated by r = ar. Therefore, the 
difference Ar(r') = r(ro + r', 6q) — r(ro, 6q) satisfies 

^oe«Ar(r') _ 

which implies Ar(r') Ri r'/arQ. Hence, 


g-Cr(ro+r',eo) _ g-C'r(ro,6»o) ^ ^-Crirofio) ^ 

aro' 


(24) 


Consider two nearby trajectories, 0(r;ro,0o) and 0(r;ro,0o + 9'). From flTSl) . 
0(r; ro, 0o + 9') — 9(r; ro, 0o) < 9' for all t > 0, therefore in view of (iii) 

|r(r, 0(r; ro, 0o))| [ar +/Sr^ cos30| a 

which implies that 

\T{ro,9o + 9') - T{ro,9o)\ < |30'5(0o)r(ro, 0o)|- 


Hence, 

|g-C(r(ro,eo+e')) _e-t^Ar-oyo)| < 
Together fl2Tl) and fl25l) imply (v). 


(25) 

QED 
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4.2.2 Proof of Theorem [2] 


As in the proof of Theorem [T|, we approximate trajectories in the vicinity of the cycle by 
superposition of local and global maps. We consider g = 0i'0202'0i ^ —)■ where 

the map 0 i is given by (jl]). The map 'ip2 in the leading order is linear: 

(rCl, Zi) = 'lp2{'V2i ^ 2 ) = (- 611^2 + B 12 Z 2 , B 21 V 2 + B22Z2)- (26) 

Because of (i), for small fi the expanding eigenvalue of ^2 depends linearly on /r, therefore 
without restriction of generality we can assume that 62 = Generically, all other eigenvalues 
and coefficients in the expressions for local and global maps do not vanish for sufficiently 
small /i and are of the order of one. We assume them to be constants independent of /i. 
From (ii), the eigenvalues satisfy ei > 0 , —Ci < 0 and —C2 < 0 . 

For small enough < 5 , in the scaled neighbourhoods B^{^2) the restriction of the system 
to the unstable manifold of ^2 in the leading order is s = /is + / 3 s^, where we have denoted 
s = W2 + iq2- As in the proof of Theorem [ 1 ], (pi,6*i) and (^2,6*2) denote polar coordinates 
in and H2"'\ respectively, such that Vi = pi cos 6^1, qi = pi sin 6^1, W2 = p2 cos 6*2 and 

q2 = p2 sin 62. We assume that the local bases near ,^1 and ^2 are chosen in such a way that the 
heteroclinic connection ^2 .^1 goes along the directions arg(6'j) = 0 for both j = 1, 2, which 

implies /3 > 0 . In the complement subspace the system is approximated by the contractions 
ii = —r2U and v = —C2V. In terms of the functions r(r, 9 ) and 'd(r, 9 ) introduced in Lemma 
[ 3 l the map 02 is 


{v2,q2) = h{P2,92) = {vo,2e sin?9(p2, ^' 2 )). 

In physical space, there are two heteroclinic trajectories from ,^1 to ^2, with positive or 
negative Wi, and three trajectories from ^2 to ^ 1 , with 6*2 = 0, 27r/3 or dvr/S. Let be 

the crossection of the heteroclinic trajectory ^1 —>■ ^2 with positive Wi. We take a modihed, 
compared to ([6]), map 'ipp. 


(P2, 92) = 01 (pi, 9 i) = 7(Api, 9 i + 0 ), where 7 = k'p^- ( 27 ) 

By choosing s, we obtain 

(p2, 6*2) = 0 i(pi, 9 i) = (Api, (— 1 )*©'), where 0 ' = 0 + 27 rs /3 satishes — 7r/3 < 0 ' < 7r/3. 

(28) 

According to Lemma [ 3 ](iv) , for small p2 and p 

sin'9{p2, 92 ), 

therefore we have 02(1^2, <?2) ~ (7?i25'2, 7?22<?2)- For small 6*1 we have sin6*1 ~ tan6*1 ^ 6*1. 
Taking into account (jl]), (| 28 |) . Lemma [ 2 ] and Lemma [ 3 ](iii), we obtain that 

g{pi,9,) ^ fc'i(piA0sinF3 30 + ^5(0))3'=iM^C'2(piA0sinF3 30 + ;,^(0))Y ( 29 ) 
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where Ci = r;o,i3 ^(3 and C 2 = (—l)^no,i3 ^(3 ^ 822 - (Here, we have denoted 

0 = |0'| and the power I in fl28l) is chosen in snch a way that in (tci, qi) = V’ 202 ' 0 i(pi 5 0) the 
value of wi is positive for small pi.) 

(a) From (1^ . the p-component of g satisfies 

Pp(pi, 0 i) > where ^3 = C'i(H/3sinV3 0)3'=iM, 

hence if 3ci < ei then for any 0 < 5 the iterates p"(pi,6'i) with initial 

0 < Pi < 5 satisfy p” (pi, 6 ^ 1 ) > h for sufficiently large n. 

(b) Assume that 3ci > ei. Since pi = 0 is a solution to 

Pp(pi,^i)-Pi = 0 (30) 

for p = 0, by the implicit function theorem, the equation fl5(l has a unique solution pi(p), 
0 < p < pi, for some pi > 0 and dri/dp = at p = 0. (In 

order to apply the implicit function theorem, we dehne the function h(pi, p) = pi — Pp(pi, 6 ^ 1 ) 
for negative pi by setting h(—pi,p) = —h.(pi,p). Note, that Pp(pi, 6 'i) does not depend 
on 6 ^ 1 .) Therefore, for small p the hxed point of the map g fl2^ can be approximated by 
{Pp,9p) = (C'i(pS'(0))3'^i/ei^ (^ 2 pS'( 0 )). This hxed point is an intersection of a periodic orbit 
with The distance from (pp. Op) to X depends on p as therefore the trajectory 

approaches X as p —)■ 0. For a given small h > 0, to hnd p{6) we approximate trajectories 

near ki and K 2 by linear (in (ni,gi) and (v 2 ,q 2 ), respectively) maps, near ^1 we use the 
approximation ([2]) and near ^2 we employ Lemma [2J We do not go into details. 

To study stability of the hxed point (pp, Op), we consider the diherence p(pi + r', 6*1 + O') — 
p(pi, 6 'i), assuming that r' and 6 ' are small and (pi, 6 'i) are close to {pp,0p). For the maps 
- 01 , '02 and 01 we have 

01 (pi + r', 0i + O') - 01 (pi, Oi) = {Ar', O'), 

02(^2 + V', q 2 + q') - 02(t^2, ^ 2 ) = (^110 + Bi 2 q', B 21 V' + ^220), (31) 

0i(wi + w', qi + q') - 01 (wi, gi) = ((no,iCi/ei)u;J^^''^"^w', (1 + (gi/wi)^)“^^0i0)- 

Recall that (see Lemma [2](v)) 

|g-C 2 (r(p 2 +r'P 2 + 0 ')) _ e-C 2 r(p 2 ,e 2 )| < g-C 2 r(p 2 , 62 ) | 0(1 + 30^(02)t(P2 , ^ 2 ) | • (32) 

PP 2 

From (l2Ti) . for small p and P 2 <C p 

I'd(p 2 + r', 62 + O') - ' 0 {p 2 , 6 ^ 2)1 ^ 'd 0 ^(p 2 , 6'2)(r'sin0^30 + ^psin”^^ 3Q//3). (33) 

By arguments similar to the ones applied in the proof of Lemma [3](iv), for small p and p 2 
the r.h.s. of fl32l) is asymptotically smaller than the r.h.s. of fl33|) . Therefore, combining fl33l) 
with (|3T1) . we obtain that for small p and pi <C p 

p(pi + r', 0i + O') - p(pi, 0i) ^ (Aip(3-i-^i)0V' + A2p3'^i0i0, AapV' + A^p^O'), 
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where the constants Aj depend on no,i, Ci, Ci, / 9 , -B12, -B22 and 0 . Since 3 ei > Ci, for small 
fi we have 

\g{pi + r\ei + e') -g{pi,ei)\ < p‘^\{r\e')\, 

and q = min(( 3 ci — ei)/( 2 ei), 1 ), which implies that the bifnrcating periodic orbit is asymp¬ 
totically stable. QED 


Remark 3 In the proof of the Theorem we have considered the map g : —)• ^ 

g = 0i'0202'01; where fji involves the symmetry 7 = tdp^ and the choice of I and s depends on 
0 and Bi 2- Hence, depending on the values of these constants, there can exist geometrically 
different periodic orbits in the vicinity of the group orbit of the heteroclinic cycle —>■ ^2 — 

with different number of symmetric copies of heteroclinic connections Ki and K2- The cycle 
shown in Fig. involves two symmetric copies of the connection ki and one connection K2. 
For a different group T (ZSO( 4 ) considered in Section\^ we present examples of various 
periodic orbits which are obtained by varying the coefficients of the respective normal form, 
see Fig. 


4.3 A numerical example 

In this snbsection we present an example of a F-epnivariant vector held possessing a hetero¬ 
clinic cycle, introdnced in proposition [H where the expanding eigenvalne of df{f2) is small. 
In agreement with Theorem O asymptotically stable periodic orbits exist in the vicinity of 
the cycle. To constrnct the nnmerical example, we start from a Lemma that determines the 
strnctnre of F-eqnivariant vector helds (the proof is left to the reader): 

Lemma 4 Any T -equivariant vector field is 


Zi = A, 
Z2 = B, 


nin2R3R4^1 


-712 -714 

/C-i /VO 




-712 ^714 


2:1 "2:0'" Z 


2 ? 


( 34 ) 


where and i?„in2n3n4 are real, 7^ 0 whenever = 0 (mod 3 ) and 

Bnin2n3n4 7^ 0 whcnevcr — n4 = l(mod 3 ). 

Keeping in fIM)) all terms of order one and two, and several terms of the third order, we 
consider the following vector held: 


Zi — aiZi 02^1 + a^zf -t- a 4 zf -t- 05^1^1 -t- a^Z 2 Z 2 + O'jzfz-i F a%Z\Z 2 Z 2 + 092:2 + aiQZ^ 
Z2 = biZ 2 + 62^1^2 + ^ 3 ^ 1^2 + & 4^2 “I" ^5^22:1^1 + &6^2^2) 


( 35 ) 


where 


we set 03 -|- 04 -|- 05 = 0 , 07 = Og = — 1 , and denote a = Oi -|- 02, a' = ai — 02- ( 36 ) 


The system fIMl) restricted into L = FixDg is 

xi = axi — xf, 
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which implies that the steady states with xi = and ^2 = —^1 exist whenever a > 0. 
The radial eigenvalnes of are — 2 a. The non-radial eigenvalnes of d/(^i) are 

—Cl = ai — 02 + 2(03 — 04) — a along the eigendirection yi , , 

Cl = + (&2 + &3)\/a + along the eigendirections X2 and 1/2 ^ 


and the eigenvalnes at df{^2) are 


02 = Oi — 02 — 2(03 — a4)y/a — a along the eigendirection yi 
—C2 = bi — (62 + h^)\fa + 65a along the eigendirections X2 and y2 

The eigenvalnes ci and —C2 are donble. 

With fl 36 p taken into acconnt, the restriction of the system 0351 ) into the plane Pi is 


Xi = axi — 2(03 + ai)y‘l — xi{x\ + yi) — XiX^ 

Pi = a'ni + 2(03 - a 4 )xipi - pi (xf + af) 


(39) 


and the one into the plane P2 is 

xi = axi + 06X2 — x\ — xixl + (09 + aiii)x\ 

X2 = biX2 + (62 + b3)xiX2 + bixl + b5x1x2 + bexl 


(40) 


In both planes the system has a form which is well-known to prodnce saddle-sink connections 
between eqnilibria and ^2 if certain conditions are fnlfilled, see e.g. p. A snfficient 
condition for existence of a heteroclinic connection ^1 —)■ ^2 in Pi is 

O3 -|- CI4 > 0 , CI3 — CI4 > 0 , CI3 — 04 -|- ((03 — ^4)^ T CK . ( 41 ) 


The expression for a condition of the existence of a connection ^2 ^1 in P2 is too bnlky, 

and therefore is not presented. 

We choose the valnes of the coefficients so that the conditions for existence of heteroclinic 
connections in Pi and P2 are satisfied, with ,^1 being a saddle in Pi and a sink in P2 and vice 
versa for ,^2- We also choose the coefficients snch that the stability condition 3 ci > ei (see 
Theorem [ 2 ]) is satisfied. 

Finally we adjnst coefficients so that the nnstable, donble eigenvalne 62 at ^2 verifies is small. 
The simnlations are performed with the following valnes of coefficients: 

a = 0 . 3 , oi = 0 . 2 , 03 = 0 . 3 , 04 = — 0 . 05 , 05 = — 0 . 25 , = 0 . 6 , 

07 = flg = — 1 , 09 = 0 . 1 , Clio — 0 . 15 ; ( 42 ) 

bi = 0 . 2 , 62 — — 0 . 1 , 63 = — 0 . 09 , 64 = — 0 . 1 , 65 = — 1 , = — 1 , 


which implies that the eigenvalnes are —ci = — 0 . 2041 , ei = 0 . 2834 , —C2 = — 0 . 4834 , 62 = 
0 . 0041 . 

In Fig. [ 2 ] we show projection of the periodic orbit and of the gronp orbit of heteroclinic 
connections, comprising the cycle, on a plane in R^, and time series of Xi and yi. The 
periodic orbit follows two heteroclinic connections ^1 ^ ^2 in Pi, and only one ^2 —t ^1 in 
P2, see remark [ 3 l 

In Section [6] we shall explore nnmerically the dynamics near the psendo-simple hetero¬ 
clinic cycle with a more involved symmetry snbgronp of SO( 4 ). 
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Figure 2: Projection of the periodic orbits (solid line) and of the heteroclinic connec¬ 
tions ^2 —^ 'Cl (dashed lines) into the plane < vi,V 2 >, where vi = (4,2,4,1.5) and 
V 2 = (2,4,—1.5,4), (a). The steady state Ci is denoted by a hollow circle and C 2 by hlled 
one. Temporal behaviour of Xi and yi, (b). 
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5 Stability of pseudo-simple heteroclinic cycles when 

r c 0(4), but r (z; so(4) 

Here we show that presence of a reflection in the group F can completely change the asymp¬ 
totic dynamics of the system ([ 1 ]) near a pseudo-simple heteroclinic cycle. As in the previous 
Section we hold on an example, which is built as an extension of the group F generated by 
transformations ([H]). Let’s introduce the reflection 

a: {zi,Z 2 ) ^ {zi,Z 2 ) (43) 

and dehne F = F U aF. This group admits the same one-dimensional and two-dimensional 
hxed point subspaces as F, therefore it also admits (similar) pseudo-simple heteroclinic cycles. 
However, it also has two types of three-dimensional hxed point subspaces: 

V = {(xi,yi,X 2 , 0 )} = Fix(cr), 

W = {(xi, 0 ,X 2 ,y 2 )} = Fix (ctk). 

Note that (i) V = Pi 0 P 2 (the invariant planes introduces in Section 14.Ih and (ii) V has 
two symmetric copies V = pV and V = p^V, while PF is a singleton. It follows from point 
(i) that if it exists, the heteroclinic cycle lies entirely in V. 


5.1 A Theorem 

Theorem 3 Consider the F -equivariant system 

X = /(x), where /(yx) = 7 /(x) for a// 7 G f (44) 

and f is a smooth map. Suppose that the system possesses equilibria ^1 and ^2 o-T^d a pseudo¬ 
simple heteroclinic cycle, introduced in proposition{J\ where the heteroclinic connection ^2 
^1 belongs to the half-plane of P 2 satisfying X 2 > 0. We assume that &^fx.^ldx\ > 0, where 
fx 2 denotes the X 2 -component of f. Let —Cj and ej be the stable and unstable eigenvalues of 
df{fj), respectively, such that —Ci and 62 have multiplicity 2 . 

Then, for sufficiently small expanding eigenvalue 62 the cycle is f.a.s. (see definition^. 

We begin with a proof of a Lemma about properties of trajectories in a generic D 3 - 
equivariant second order dynamical system ffTTj) . where we assume that a > 0 and > 0 . 
Recall, that in polar coordinates the system takes the form ffT^ . the system is considered in 
the sector 0 < 6 ^ < vr/S, by r(ro, 6 * 0 ) we have denoted the time it takes the trajectory of the 
system flT^ starting at (vq, Gq) to reach r = 1 and by 'i9(ro, do) the value of 6 ^ at r = 1 . 


Lemma 5 7/0 < 6*0 < vr /6 then i9(ro,do) satisfies 


sm3'd{ro,9o) < 


ro(d cos 300 0 a 
(3 cos 300 + 


3 

sin 300. 


(45) 
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Proof: From (fT^ . 

— (r® sin^ 36 ) = 6ar® sin^ 36 , 
at 

which implies that 


r®(t) sin^ 36 {t) = e®“V® sin^ 300 - 
Setting t = r(ro, 6 q) we obtain that 

sin 30 (r) = sin 30 o- ( 46 ) 

As well, flT^ implies that 0 < 0 and that for 0 i < 62 we have f(r, 0 i) > r(r, 62). Hence, the 
trajectory (r(t),0(t)) starting at r(0 ) = tq and 0(0) = 0o satisfies 

r(r(t),0(t)) > f(r(t),0o) = ar + /dcos 30 or^ for all t > 0. 

Therefore, by the same argnments as we employed in the proof of Lemma [ 3 ](i,ii) 

g«r(ro,eo) < rp/j cos 300 + a 
ro{(3 cos 300 + a) 

Combining fl 46 ]l and flT 7 )l we obtain statement of the Lemma. 


(47) 

QED 


Proof of Theorem [3} We aim at finding conditions snch that the first return map g = 
'0i</>i'0202 : —)■ is a contraction and hence converges to 0 as the number of iterates 

tends to infinity. (Note, that here we consider the map while in Theorems [ 1 ] 

and [ 2 ] it was —)■ The local and global maps are defined the same way, as it 

was in the proofs of Theorems [ 1 ] and [ 2 l Because of the refiection a, the constants 0 , B12 
and B21 in ([ 6 ]) and fl 26 |l vanish, hence global maps take the form 

(P2, 62) = -01(01, di) = {Api, 0 i), (wi, qi) = 'Ip2{v2, q2) = {B11V2, B22q2)- ( 48 ) 

Recall, that ( 0 i, 0 i) and (02,02) denote polar coordinates in and H2"'\ respectively, 

such that Vi = pi cos 61, qi = pi sin 61, W2 = p2 cos 62 and q2 = p2 sin 62- The map 0 i is given 
by dl]). The map 02 is 

( 112 ,^ 2 ) = 02(02,02) = (no,2e“''2^^^=’®^\tan'0(p2,02))- 
The composition g = 0i0i0202 therefore is 


9{f’2,02)^ ('l?ii%2e ,arcsin(i?22tan,?(p2,^2 


( 49 ) 


The condition d‘^fx2/dx2 > 0 implies that the restriction of the system into the 
unstable manifold of ^2, and which is of the form (fTTD . satisfies 0 > 0. The condition a > 0 
is evidently satisfied because of the instability of ^2, hence we can apply Lemmas [ 3 ] and [ 5 ] to 
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estimate t{p 2 , 62 ) and '^(^ 2 ,^ 2 )- Denote by Pp and ge the first and second components of g 
in fl49p . From Lemmas [3]^i,ii) and Owe can write the estimates 


ciC2/eie2 


ge{'f' 2 ,d 2 ) < arcsin I i? 22 tan I l/3arcsin 


f p2/3cos36>2+e2 

y ^cos36»2+e2 



(50) 


where C = Avq^i{BiiVq^ 2 Y^^^^■ We show that gp and gg are contractions under the conditions 
stated in the Theorem. 

The function gp{p 2 , O 2 ) is a smooth function of p 2 for all p 2 > 0. Its derivative is positive 
and has the expression 


9'p{p2) 


hC 


62 (/^ + 62 )^ 
(P2/3 + 62)^*+^ 






where we have denoted h = The condition h > 1 implies that lim q'(p) = 0. Therefore, 

6162 ^ ’ 

there exists p > Q such that g'p{p 2 ) < 1 and Pp is a contraction when p 2 < p. 

To prove that gg is a contraction let us assume that 62 is chosen small enough to allow for the 

approximation tan302 ~ sin302 ~ 302 (this is not necessary but simplihes the calculations). 

Then 

MO.) < ~ 

y p cos 302 + 62 


02 . 


If we chose 62 and p 2 small enough (/9 is a hxed positive constant), then the factor in front 
of 02 is smaller than 1. Therefore, if p 2 converges to 0 by iterations of g, the same is true for 
02 . Combining the conditions for Pp and gg we get the result. QED 


Remark that the convergence to 0 of 62 is slow compared to that of p 2 - 

Remark 4 Suppose that, similarly to Theorem\^ we consider a T-equivariant system which 
depends on a parameter p, possesses a heteroclinic cycle for 0 < p < po and 62 = 0 for 
p = 0. Then by Theorem there exists > 0 such that the heteroclinic cycle is f.a.s. 
for 0 < p < pi. Moreover, the condition d‘^fx^ldx\ > 0 (which is equivalent to (3 > t)) 
is always satisfied in such a system, by the reasons given in the proof of Theorem [B Note 
however, that closeness to the bifurcation point, and even smallness of 02 , are not necessary 
for fragmentary asymptotic stability of the cycle. 


5.2 A numerical illustration of Theorem [3] 

The equivariant structure of T-equivariant vector helds is easily deduced from that of T 
equivariant vector helds. 

Lemma 6 Any T-equivariant vector field has the form fl34p (see Lemma\^ with the following 
additional conditions on the coefficients. 712113114 ^111112114113 nnd .R 111112113114 -^112111113114 ■ 
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In this subsection we consider the same system flM|) as we considered in subsection 14.31 
Because of the presence of the symmetry a (see Lemma E]), coefficients of this system satisfy 
Og = Oio and 62 = & 3 - As in subsection 14.31 we denote a = Oi + 02 , a' = Oi — 02 and assume 
03 + 04 + 05 = 0 and oy = os = — 1. Therefore, the system is 

Zi = ctiZi + 02^1 + + a^zf + a^ZiZi + 05^2^2 — zfzi — T ^9(^2 T /ci \ 

Z2 = &1^2 + ^2 (^1^2 + Z1Z2) + 64^1 + 55^2^1^! + ^ 6 ^ 2^25 

The conditions for the existence of the heteroclinic cycle in the invariant space V are 
similar to those given in Section 14.31 It is however interesting to analyse the dynamics in 
the invariant space W. Consider the restriction of system fl5TD into the subspace W: 

xi = axi + 03(0:2 + yl) — xf — xi(x2 + y^) + 2ag(x2 — 3x21/2) 

±2 = biX 2 + 2 b 2 XiX 2 + ^ 4 ( 3:2 - yl) + b^X 2 x\ + beX 2 {xl + j/|) (52) 

m = biy2 + 2 b2Xiy2 - 2643 : 21/2 + 651 / 23 ^? + bey2{xl + //|). 

In the limit 64 = 0 any plane containing the axis L is flow-invariant by (|5^ . In fact, any 
such plane is a copy of P 2 by some rotation around L. Therefore, whenever a saddle-sink 
connection exists in P 2 , it generates a two-dimensional manifold of saddle-sink connections 
in W (an ’’ellipsoid” of connections). Small T-equivariant perturbations of the system, in 
particular switching on 64 7 ^ 0, do not destroy this manifold. If in addition the conditions 
an are satisfied, then a continuum of robust heteroclinic cycles connecting and ^2 exists 
for the T-equivariant vector field. 

We take the same values of the coefficients as in fl42D . except that Og = Oio = 0.1 and 
62 = 63 = —0.6 so that the system is T-equivariant. 

Note, that when the symmetry a is absent, the heteroclinic cycle is completely unstable 
but a nearby asymptotically stable periodic orbit is observed exists, in accordance with the 
results of previous Sections. When the symmetry a is present, the trajectory starting from 
some point in a neighborhood of ^2 asymptotically converges to the heteroclinic cycle (or to 
its symmetric copy) as f —)■ 00 , see Fig. |3l The convergence in the (xi,X 2 ) plane is slower 
than in the (xi,//i) plane, because the convergence 6 * —?■ 0 is slower than p —)■ 0 . 

6 Another example with a group P C *5'0(4) 

The groups which have been considered in Sections 0] and [5] had the advantage of allow¬ 
ing for an elementary algebraic description, which gives an opportunity to study in detail 
the structure on invariant subspaces and related dynamics. In this Section we aim at pro¬ 
viding an example of an admissible group with pseudo-simple heteroclinic cycles, with a 
non-elementary group structure. In principle, all such groups can be listed as we did for 
admissible groups with simple heteroclinic cycles in [T3] . 

We then build an equivariant system for this group and study numerically its asymptotic 
behavior near the pseudo-simple heteroclinic cycle. We begin with a description of the 
bi-quaternionic presentation of finite subgroups of S'0(4) together with some geometrical 
properties of these actions, which we employ to identify invariant subspaces. 
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Figure 3: Projection of a trajectory of the system (l5T]) converging to a heteroclinic cycle (solid 
line) and of heteroclinic connections (dotted, dashed and dashed-dotted lines, different types 
of lines label different group orbits of heteroclinic connections) into the plane < vi, V 2 >, 
where vi = (4, —2,4,1.5) and V 2 = (2,4, —1.5,4), (a). Projection of this trajectory into the 
planes [xi^yi) (b) and {xi,X 2 ) (c). 
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6.1 Subgroups of 50(4): presentation and notation 

Here we briefly introduce the bi-quaternionic presentation and notation for finite subgroups of 
S'0(4), which has been used in [13] to determine all admissible groups for simple heteroclinic 
cycles in See [3] for more details. 

A real quaternion is a set of four real numbers, q = {qi,q 2 ,(l 3 ,Qi)- Multiplication of 
quaternions is defined by the relation 

qw = (giwi - q 2 W 2 - qsWs - q^w^, qiW 2 + q 2 Wi + q^Wi - q^w^, . . 

qiwz - q2W4 + qswi - ^ 4 ^ 2 , qiW4 + ^ 2^3 - q3W2 + q4Wi). 

q = (gi, —q2, —qs, —^4) is the conjugate of q, and qq = g? + gl + + ?! = IhP is the square 

of the norm of q. Hence q is also the inverse q“^ of a unit quaternion q. We denote by Q 
the multiplicative group of unit quaternions; obviously, its unity element is (1,0, 0,0). 

Due to existence of a 2-to-l homomorphism of Q onto SO(3) (see [3]), finite subgroups 
of Q are named after respective subgroups of SO (3). They are: 

= ©”“q (cos2r7r/n, 0, 0, sin2r7r/n) , D„ = Z 2 „, © ©^4o^(0, cosrvr/n, sinrTr/n, 0) 

V = ((±1,0,0,0)) , T = V© (±i, ±i, ±|, ±i) 

O = T© y4((±l, ±1,0,0)) , I = T©4((±r,±l,±r-S0)) 

(54) 

where r = (\/5 ± l)/2. Double parenthesis denote all possible permutations of quantities 
within the parenthesis and for I only even permutations of (±r, ±1, ±r“^, 0) are elements of 
the group. 

By regarding the four numbers (gi, g 2 , gs, g 4 ) as Euclidean coordinates of a point in R"^, 
any pair of unit quaternions ( 1 , r) can be related to the transformation q —lqr~^, which is an 
element of the group SO(4). The respective mapping $ : Q x Q —SO(4) is a homomorphism 
onto, whose kernel consists of two elements, (1,1) and (—1, —1). Therefore, a hnite subgroup 
of SO(4) is a subgroup of a product of two hnite subgroups of Q. There is however an 
additional subtlety. Let T be a hnite subgroup of SO(4) and Q = <l)“^(r) = {(b,rj), 1 < 
j < J}- Denote by L and R the hnite subgroups of Q generated by \j and r^, 1 < j < J, 
respectively. To any element 1 G L there are several corresponding elements r*, such that 
(1, Fj) G Q, and similarly for any r G R. This establishes a correspondence between L and 
R. The subgroups of L and R corresponding to the unit elements in R and L are denoted 
by Ijk and Hk, respectively, and the group ^ by (L | L^; R | R_ft:). 

Lemma 7 (see proof in Let Ni and N 2 he two planes in R"^ and pj, j = 1,2, be the 
elements of SO(4) which act on Nj as identity, and on as —I, and ^~^pj = (1^, r^), where 
$ is the homomorphism defined above. Denote by (lil 2 )i and (rir 2 )i the first components of 
the respective quaternion products. The planes Ni and N 2 intersect if and only if (lil 2 )i = 
(rir 2 )i = cos a and a is the angle between the planes. 
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Lemma 8 Consider g E SO(4), ^ = ((cosa,sinav); (cos/9,sin/9w)). 

Then dim Fix < g >= 2 if and only if cos a = cos/3. 

Lemma 9 Consider g^s eS 0(4), where ^~^g = ((cos a, sin av); (cos a, sin aw)) and 
$-is = (( 0 ,v); ( 0 ,w)). 

Then Fix < g >= Fix < s >. 


6.2 The group F = (D3 | Z2; O | V) 

The group (D 3 |Z 2 ; 0 |V) was proposed in [TB] as an example of a subgroup of 0(4) ad¬ 
mitting pseudo-simple heteroclinic cycles. It is the (unique) four dimensional irreducible 
representation of the group GL{2, 3) (2x2 invertible matrices over the held Z 3 ). This group 
is generated by the elements cxi (order 8) and (J 2 (order 2 ) below: 


0-1 


02 \ _/^ 20 \ 
2 2 ) ^ yo 1 ) 


In quaternionic form its elements are: 


(( 1 , 0 , 0 , 0 ) ; 
((cos( 7 r/ 3 ), 0, 0, sin( 7 r/ 3 )) ; 
((cos(27r/3), 0, 0, sin(27r/3)) ; 

(( 0 , 1 , 0 , 0 ) ; 
((0, cos( 7 r/ 3 ), sin( 7 r/ 3 ), 0) ; 
((0, cos(27r/3), sin(27r/3), 0) ; 


V), 

(l,l,l,l)/2 V), 
(-l,l,l,l)/2 V), 
(l,l,0,0)/x^ V), 
(l,0,0,l)/^/2 V), 
(l,0,l,0)/^/2 V). 


Lemmas [H] and [9] imply that the group has two isotropy types of subgroups S satisfying 
dim Fix S = 2. They are generated either by an order three element, 

e(s, r) = ((cos( 7 r/ 3 ), 0, 0, sin( 7 r/ 3 )); (1, (-1)", (-1)’’, (-l)"+’’)/2), 

where r, s = 0,1, (the proof that dim Fix < e >= 2 follows from Lemma [ 8 ]) or by an order 
two element, a plane rehection 

K{q, n, t) = (- 1 )''(( 0 , cos(n 7 r/ 3 ), sin(n 7 r/ 3 ), 0 );p’"(( 0 , 0 , (- 1 )*, l)/V 2 )), 


where n = 0,1, 2, g, t = 0,1 and p is the permutation p(a, b, c, d) = (a, c, d, b). The isotropy 
planes can be labelled as follows: 


Pi(s,r) = FixSi(s,r), where Si(s,r) =< e(s,r) >, 

P 2 {q,n,t) = FixS 2 (g,n,f), where T, 2 {q,n,t) =< > • 

Each of these plane contains exactly one copy of each of two types of (not conjugate) sym¬ 
metry axes, the isotropy groups of which are isomorphic to D 3 . They are: 


Li(s, r) is the intersection of Pi(s, r), P 2 ( 0 , 0, s -|- 1 ), P 2 (s -|- 1 , 1 , r) and P 2 (r’, 2 , r -|- 1 ), 
L 2 {s,r) is the intersection of Pi(s, r), P 2 (l, 0 , s-|- 1 ), P 2 (s, 1 , r) and P 2 (r-|- 1 , 2 , r-|- 1 ). 
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(Intersections of Pi and P 2 or of two P 2 S follows from Lemma [71) Since iV(Ei)/Ei = D 2 , 
the axes Li(s,r) and L2{s,r) intersect orthogonally. 

Now, if 

(i) there exist two steady states ^i E Li and ^2 ^ L 2 ; 

(ii) .^1 is a saddle and ^2 is a sink in Pi, moreover a saddle-sink heteroclinic orbit ki exists 
between ,^1 and ^2 in Pi; 

(iii) ^1 is a sink and ^2 is a saddle in P 2 , moreover a saddle-sink heteroclinic orbit K 2 exists 
between ^1 and ^2 in Pi, 

then there exists a psendo-simple robust heteroclinic cycle ^2 —t ^1 ^ ^ 2 - If in addition the 
unstable manifold of ^2 in P 2 is included in the stable manifold of (hh“(^ 2 )nPi) C hL^(^i), 
then the system possesses a heteroclinic network comprised of two distinct (i.e. not related by 
any symmetry) connections ^2 ^1 and one connection ^1 ^2- With minor modihcations 

of Lemma 5 in [13], it can be proven that the conditions (i)-(iii) above are satished for 
an open set of L-equivariant vector helds. In the following subsection we build an explicit 
system with these properties. 

6.3 Equivariant structure 

Here we derive a third order normal form commuting with the action of the group L = 
(D 3 I Z 2 ; O I V) introduced in subsection 16.21 For convenience, we write the quaternion q = 
(gi, q2, qs, ^ 4 ) as a pair of complex numbers, 

u = (mi, U2 ), where Mi = gi -|- ig 2 and M 2 = <?3 + (55) 

The operation of multiplication fl53p takes the form 

hu = {hiUi - h2U2, hiU2 h2Ui) (56) 

and u = (hi, —M2) is the conjugate of u. For quaternions presented as ( 1551 ) and points in 
as z = (zi, Z2), the action ( 1 ; r) : z ^ Izr^^ of (some) elements of F is 


(l;r) 

z —)■ Izr ^ 

((0,l,0,0);(l,l,0,0)/^/2) 
«(l,0,l) = ((0,l,0,0);(0,0,l,-l)/x/2) 
«(0,0,l) = ((0,l,0,0);(0,0,-l,l)/x^) 
6(0,0) = ((l,0,0,^/3)/2;(l,l,l,l)/2) 

{zi,Z2) -E (eWSi,e3WS2) 

(^1,^2) (e^W4^2,e“^W4^^) 

{zi,Z2) -E -(e3W4^2,e-3W4^^) 

(zi, Z 2 ) -E (e-W4^^ + e-W4^2 - i\/3(e-W4^^ + 

-eW4;2^ + eW4^2 + i\/3(eW4^^ + 0^74^^). 


( 57 ) 
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Using the computer algebra program GAP, we derive from fl57|) that the cubic equivariant 
terms are of the form 


Zl 

Z 2 


b{zizi + Z2Z2) 


Zl 

Z 2 


+ C01 + (i02 + 603 


where 6, c, d, e are real and the maps 0j are cubic expressions of zi, Zi, Z 2 , Z 2 '- 


f \Piz\z 2 - ^zfzi + ^ 1 ^ 2-22 - J^1^2 \ 

\ -\/3ZiZ2 - ^Z2Z2 + iZiZiZ2 - fzfz2 J 


a, _ f ^2 ~ ^^ 1^2 + 2i^/3zlZ2 \ 

~ ^ V -zf + 82 : 1^2 + 2i\/3zfz2 J 

^ _ f -^y/^zfzi - \/3 ziZ 2 + iz2 + 3izlz2 \ 
\ —3\/32;|z2 — '/Szlz 2 — izf — Siziz^ J 

In another arrangement of the cubic monomials, the equations read: 


Zi — (h" A bi^^iZi + Z2Z2))Zi + Azfzi + iC2:iZ2^2 “1“ Bz2 + Czfz2 + DZ 1 Z 2 

Z2 = {fi + b{ziZi+ Z2Z2))Z2 + Az2Z2 + icZiZiZ2- Bzf - CZ1Z2+ DzIz2 


(58) 


where 

A = —SVSe ~ B = d + ie, C = —3d + y/Sc + 3ie, D = \/3(—e + i{2d — ^^c)). (59) 

We set /r = 1 in order to make the origin unstable. 

From (1571) we obtain that the invariant planes are: 

P 2 (l, 0,1) = Fix (((0,1, 0, 0); (0, 0,1, -1)/V^)) = (e^^/^ra, w) 

P2(0, 0,1) = Fix (((0,1, 0, 0); (0, 0, -1,1)/V^)) = (-e^^'/^w;, w) 

Pi(0, 0) = Fix (((1, 0, 0, \/3)/2; (1,1,1, l)/2)) = {y/ 2 e-^^^Byj _ 


and the invariant axes are given by 

Pi(0,0) = Pi(0,0)nP2(0,0,l) = (-e^"/Vie“brie“i), where e^^^ 

and 

L2(0,0) = Pi(0,0) nP2(l,0,l) = (e2"/V2e"Ur2e"2), where e^"^ 


The system (1581) restricted onto the axes Pi is 


ri = ri(p + rl{2b + 2\/2c — 


8 V 2 

Tf 



-l + \/2i 




1 + \/2i i/4 


25 





( 60 ) 


and onto L 2 


Therefore, the conditions for existence of the steady states G Li and .^2 £ L 2 are 




respectively. Whenever the steady states exist, the eigenvalnes of df{^j) in the radial direc¬ 
tions are —2/i. 

By considering the restriction of the system into the snbspace Fix (—((0,1, 0, 0); (0, 0,1, —l)/\/2)) 
we obtain that the eigenvalnes of d/(^i) and d/(^ 2 ), associated with two-dimensional eigenspaces, 
are 



By considering the restriction of the system into the snbspace hxed by ((1, 0, 0, \/3)/2; (l,l,l,l)/2) 
we derive the expressions for the remaining eigenvalnes: 



( 62 ) 


Necessary conditions for the existence of heteroclinic cycle —>■ ^2 —^ ^1 discnssed in snb- 

section 16.21 are that 


Af > 0, \f < 0, A 


.f* < 0 and A“^* > 0. 


(63) 


6.4 The numerical simulations 

The general third order T-eqnivariant system is given by dSHD-dSlD- From the nec¬ 

essary conditions for existence of the heteroclinic cycle expressed in terms of the coefficients 
6, c, d and e are: 



We set 


where — 1 < hi < 1 and ^2 is positive. 
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In Section m we have proved that in a F-equivariant system, for the considered group F = D 3 , 
an asymptotically stable periodic orbit bifurcates from a pseudo-simple heteroclinic cycle 
when the double expanding eigenvalue vanishes. It can be easily shown that the expressions 
for local and global maps near a heteroclinic cycle in a (D 3 | Z 2 ; O | V)-equivariant system 
are identical to the ones given in the proof of Theorem [2], with 7 G (D 3 | Z 2 ; O | V) in ([27]). 
Therefore, the proof holds true for F = (©3 | Z 2 ; O | V) that we consider here. The expanding 
double eigenvalue is = 2^h2rl/3, hence we consider small values of h 2 . 

The computations were done for several values of hi, 0.5 < hi < 1. For small enough ^2 
the attractor is a periodic orbit near the group orbit of heteroclinic connections Ki and k, 2 , 

X = (U^er7«i) 

Three instances of these orbits are shown in Figs. Il|and[5l Note, that by varying hi we get 
periodic orbits with different 7 ’s in (127|) (see remark [3]). 

With the increase of h 2 the behaviour ceases to be time-periodic, however the chaotic 
trajectory stays close to the heteroclinic cycle. The maximal distance of the points of the 
trajectory from the heteroclinic cycle increases with h 2 (see Figs. E]and[7|). 

7 Conclusion 

Heteroclinic cycles in have been extensively studied in the last three decades, however only 
recently ra existence of pseudo-simple heteroclinic cycles has been noticed. The present 
paper is the hrst one where properties of these cycles are investigated. Our main result is 
that pseudo-simple heteroclinic cycles in are completely unstable under generic conditions 
when the symmetry group F of the system is contained in SO(4). We further analysed 
the asymptotic behavior near a pseudo-simple cycle with a specihc group F such that the 
heteroclinic cycle connects precisely two different equilibria. We proved that: (i) when the 
unstable double eigenvalue at one equilibrium is close to 0 , a nearby stable periodic orbit can 
exist under generic conditions, (ii) when F is extended to a group F which is not a subgroup 
of S'0(4), then the pseudo-simple heteroclinic cycles are fragmentarily asymptotically stable. 
These properties have been illustrated numerically. Finally, in the last part of this study 
we have considered a more complex subgroup of SO(4), for which the proof of existence of 
periodic orbits near the heteroclinic cycles is still valid but which possess a richer structure. 
Numerical simulations have shown that periodic orbits can follow different connections along 
the group of heteroclinic cycle and that non-periodic attractors in the vicinity of the cycle 
can also exist. 

We expect that the results of Section HI can be generalized to other subgroups F C SO (4), 
at least when the unstable two-dimensional manifolds are invariant by the action of © 3 . We 
also intend to look at the case when the two-dimensional unstable manifolds are invariant 
under the action of ©^ with k > 3, the considered example being readily modihed to F = D^. 
Arguments of Theorem [3] are likely to hold true for other subgroups of O (4), which are not 
in SO(4), admitting pseudo-simple heteroclinic cycles. In light of our numerical observations 
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(a) 




Figure 4: Projection of the periodic orbits (solid line) and of the heteroclinic connections 
(dashed and dotted lines) into the plane < vi,V 2 >, where vi = (1,2,1,1.8) and V 2 = 
(2, —1,1.8, —1), for /i 2 = 0.001 and hi = 0.7 (a), hi = 0.8 (b) and hi = 0.92 (c). Steady 
states ^1 are denoted by hollow circles and ^2 by hlled ones. 
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Figure 6: Projection of the chaotic trajectories (solid line) and of the heteroclinic connections 
(dashed and dotted lines) into the plane < vi,V 2 >, where vi = (1,2,1,1.8) and V 2 = 
(2,-l,1.8,-l), for hi = 0.92 and /i 2 = 0.0015 (a), /i 2 = 0.002 (b) and ^2 = 0.0028 (c). 
Steady states are denoted by hollow circles and ^2 by filled ones. 
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in Section 6 , we think it will be of interest to investigate further the transition to complex 
dynamics near a pseudo-simple cycle with symmetry group in SO(4), for groups studied in 
this paper, of for different ones. 

The dehnition of pseudo-simple heteroclinic cycles can be generalized to with n > 4 by 
requiring the unstable eigenvalue at one of the equilibria to have dimension of the associated 
eigenspace to be larger than one. Evidently, such cycles are not asymptotically stable. The 
question whether they can be f.a.s. for T C SO(? 7 ,) is yet another open problem. 


References 

[1] P. Chossat. The bifurcation of heteroclinic cycles in systems of hydrodynamical type. 
Journal on Continuous, Discrete and Impulsive System, 8a, 4, 575-590 (2001). 

[2] P. Chossat and R. Lauterbach. Methods in Equivariant Bifurcations and Dynamical 
Systems. World Scientific Publishing Company, 2000. 

[3] P. Du Val. Homographies, Quaternions and Rotations. OUP: Oxford, 1964. 

[4] G. Faye, P. Chossat. Bifurcation diagrams and heteroclinic networks of octagonal H- 
planforms, J. of Nonlinear Science, 22, 1, 277-326 (2012). 

[5] M. Field. Lectures on Bifurcations, Dynamics, and Symmetry. Pitman Research Notes 
in Math. Series 356, Longman, 1996 

[ 6 ] M. Koenig. Linearization of vector helds on the orbit space of the action of a compact 
Lie group. Proc. Camh. Phil. Soc. 121, 401-424 (1997). 

[7] M. Krupa. Robust heteroclinic cycles. J. Nonlinear Science, 7, 129-176 (1997). 

[ 8 ] M. Krupa and 1. Melbourne. Asymptotic stability of heteroclinic cycles in systems with 
symmetry. Ergodic Theory Dyn. Syst. 15, 121-148 (1995). 

[9] M. Krupa and 1. Melbourne. Asymptotic stability of heteroclinic cycles in systems with 
symmetry. If. Proc. Roy. Soc. Edinburgh 134A, 1177-1197 (2004). 

[10] S. Lang. Algebra, third edition. Addison-Wesley, 1993. 

[11] O.M. Podvigina. Stability and bifurcations of heteroclinic cycles of type Z. Nonlinearity 
25, 1887-1917, arXiv: 1108.4204 [nlin.CD] (2012). 

[12] O.M. Podvigina. Classihcation and stability of simple homoclinic cycles in Nonlin- 
eanty 26, 1501-1528, arXiv:1207.6609 [nlin.CD] (2013). 

[13] O.M. Podvigina and P. Chossat. Simple heteroclinic cycles in Nonlinearity 28, 901- 
926, arXiv: 1310:0298 [nlin.CD] (2015). 


32 


[14] O.M. Podvigina and P. Ashwin. On local attraction properties and a stability index for 
lieteroclinic connections. Nonlinearity 24, 887-929, arXiv:1008.3063 [nlin.CD] (2011). 

[15] G. W. Schwarz. Lifting smooth homotopies of orbit spaces. Publ. Math. I.H.E.S. 51, 
37-135 (1980). 


33 


